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We present a model based on the SU(3)c ® SU(3)l ® U(l)x gauge symmetry having an extra 
T 70 Z 2 ®^ 3 ®Zi 4 flavor group, which successfully describes the observed SM fermion mass and mixing 
pattern. In this framework, the light active neutrino masses arise via double seesaw mechanism and 
the observed charged fermion mass and quark mixing hierarchy is a consequence of the Z 2 ®Zs® Z \4 
symmetry breaking at very high energy. In our minimal TV flavor 331 model, the spectrum of 
neutrinos includes very light active neutrinos as well as heavy and very heavy sterile neutrinos. The 
model has in total 16 effective free parameters, which are fitted to reproduce the experimental values 
of the 18 physical observables in the quark and lepton sectors. The obtained physical observables for 
both quark and lepton sectors are compatible with their experimental values. The model predicts 
the effective Majorana neutrino mass parameter of neutrinoless double beta decay to be mpp = 
3 and 40 meV for the normal and the inverted neutrino spectrum, respectively. Furthermore, our 
model features a vanishing leptonic Dirac CP violating phase. 


I. INTRODUCTION 


In spite of the experimental confirmation of the big accomplishments of the Standard Model (SM) in describing 
electroweak phenomena, given by the discovery of the ~ 126 GeV Higgs boson by ATLAS and CMS collaborations 
at the CERN Large Hadron Collider (LHC) [JQ, there are several issues that the SM is unable to address. Some of 
these issues are the Dark Matter problem, the fermion mass and mixing hierarchy and the neutrino oscillations [5j, ib}. 
Furthermore, there are problems with the matter-antimatter asymmetry related with new phases responsible for CP 
violation. Moreover, the SM does not explain the tiny value of the neutron dipole moment. Because of these reasons 
it is necesary to consider an extension of the Standard Model. In particular, the observed the quark mass and mixing 
pattern, which cannot be explained from first principles in the context of the Standard Model, is a clear indication 
of physics beyond the Standard Model. On the other hand, the discovery of the Higgs boson opens the possibility 
to unravel the Electroweak Symmetry Breaking (EWSB) mechanism and motivates to study extensions of the SM 
having extra scalar particles that could provide an explanation for the existence of Dark Matter [j]. 

The lack of predictivity of the Standard Model Yukawa sector, motivates to consider extensions of the Standard 
Model aimed to address its flavor puzzle. Discrete flavor symmetries are important because they generate fermion 
textures useful to explain the three generation flavor structure, for recent reviews see Refs. [8H10|. These discrete 
flavour symmetries have been employed in extensions of the Standard Model with the aim to study the fermion mass 
and mixing hierarchy in order to address the flavor puzzle of the SM. Discrete flavor symmetries can arise from the 
underlying theory, e.g., string theory or compactification via orbifolding. In particular, from heterotic orbifold models, 
one can generate the D 4 and A(54) flavor symmetries 0 - Furthermore, magnetized/intersecting D-brane models can 
generate the A (27) flavor symmetry 0 - Discrete symmetries may link the low energy physics and the underlying 
theory. 


Furthermore, another unaswered issue in particle physics is the existence of three generations of fermions at low 
energies. The mixing patterns of leptons and quarks are significantly different; while in the quark sector, the mixing 
angles are small, in the lepton sector two of the mixing angles are large and one is small. Models having SU(3)c < 8 > 
SU(3 )as a gauge symmetry, are vectorlike with three fermion generations and thus do not contain anomalies 
12 416|. Defining the electric charge as the linear combination of the X 3 and Tg SU(3)l generators, we have that it 
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is a free parameter, which does not depend on the anomalies (/3). The charge of the exotic particles is defined by 
setting a value for the /? parameter. Setting /3 = — -^=, implies that the third component of the weak lepton triplet is a 
neutral field v ^ allowing to build the Dirac Yukawa term with the usual field of the weak doublet. By adding very 
heavy sterile neutrinos iV ^’ 2 ' 3 in the model, light neutrino masses can be generated via double seesaw mechanism. The 
331 models with j3 = — ^ provide an alternative neutrino mass generation mechanism and include in their neutrino 
spectrum light active sub-eV scale neutrinos as well as sterile neutrinos which could be dark matter candidates if they 
are light enough or candidates for detection at the LHC, if they have TeV scale masses. Having TeV scale sterile 
neutrinos in its neutrino spectrum, makes the 331 models very important since if these sterile neutrinos are detected 
at the LHC, these models can shed light in the understanding of the electroweak symmetry breaking mechanism. 

Neutrino oscillation experiments [§, TtHU show that there are at most one massless active neutrino and that 
the different neutrino flavors mix. Neutrino oscillations experiments do not determine neither the absolute value of 
the neutrino masses nor the Major ana or Dirac feature of the neutrino. Nevertheless neutrino mass bounds can be 
obtained from tritio beta decay 22 |, double beta decay [23| and cosmology 24 1 . 

The global fits of the available data from the Daya Bay H3, T2K |l 8 j], MINOS pj|, Double CHOOZ [2(| and 
RENO [21] neutrino oscillation experiments, constrain the neutrino mass squared splittings and mixing parameters 
[25l |. The current experimental data on neutrino mixing parameters suggests a violation of the tribimaximal symmetry 
described by the Tribimaximal Mixing (TBM) matrix, whose predicted mixing angles satisfy (sin 2 0 i 2 ) TBM = 
(sin 2 6 * 23 )„ D .. = A, and (sin 2 = 0. To generate nearly tribimaximal leptonic mixing angles consistent 

' ' 1 hi IVl Z , ' ' 1 jd 1V1 1 1 I 1 

with the experimental data, discrete symmetry groups [26H35| are implemented in extensions of the Standard Model. 
Another approach to address the flavor puzzle consists in postulating fermion mass textures (see Ref [3(| for some 
works considering textures). Moreover, models based on extended symmetries in the context of Multi-Higgs sectors, 
Grand Unification, Extradimensions and Superstrings have been explored [§, [37ti4fij to provide an explanation for 
the observed fermion mass and mixing pattern. Furthermore, in the framework of minimal 331 models, the discrete 
groups require of an extra high energy scale, larger than the scale of breaking of the SU(3)c ® SU(3 )l ® U{l)x 
symmetry, and in some cases, new scalar fields need to be introduced at the very high discrete symmetry breaking 
scale, with the aim to fulfill the irreducible representations of these discrete groups that allow to get viable fermion 
textures at low energies, after the gauge and discrete symmetries are spontaneosly broken. 

In this paper we formulate an extension of the minimal SU(3)c x SU(3)l x U(1)x model with /3 = —- 7 =, where 
an extra T 7 ® Z 2 ® Z 3 ® Z 14 discrete group (see Ref [35| for studies about the X 7 flavor group) extends the symmetry 
of the model and very heavy extra scalar fields are added with the aim to generate viable and predictive textures for 
the fermion sector that successfully describe the fermion mass and mixing pattern. The obtained physical observables 
in the quark and lepton sector are consistent with the experimental data. Our model at low energies reduces to the 
minimal 331 model with /3 = — ^=. 

The content of this paper goes as follows. In Sec. [TT] we describe the proposed model. Sec. eh is devoted to the 
discussion of lepton masses and mixings. In Sec. El we present our results in terms of quark masses and mixing, 
which is followed by a numerical analysis. Conclusions are stated in Sec. |Vj Appendix |A] provides a brief description 
of the T 7 discrete group. 


II. THE MODEL 

A. Particle content 


We consider a SU(3)c < 8 > SU(3)l ® U( l)x ® T 7 ® Z 2 ® Z 3 ® Z 14 model where the full symmetry Q is spontaneously 
broken in three steps as follows: 


g = SU(3) C <8> SU (3) l <g> U (l) x Z 14 ^+ 

SU(3) C ® SU (3) L (8) U (l) x ^SU{3) C <8> SU (2) L ® U (1V 
SU(3) c ®U( 1) q , 

where the hierarchy v v , v p <C v x <C A m j among the symmetry breaking scales is fullfilled. 
The electric charge in our 331 model is defined as: 


Q = t 3 


T3 TS + XI ’ 


(1) 


( 2 ) 
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where X 3 and Tg are the SU(3)l diagonal generators and I is the 3x3 identity matrix. 

Two families of quarks are grouped in a 3* irreducible representations (irreps), as required from the SU(3)l anomaly 
cancellation. Furthermore, from the quark colors, it follows that the number of 3* irreducible representations is six. 
The other family of quarks is grouped in a 3 irreducible representation. Moreover, there are six 3 irreps taking into 
account the three families of leptons. Consequently, the SU(3)l representations are vector like and do not contain 
anomalies. The quantum numbers for the fermion families are assigned in such a way that the combination of the 
U(l)x representations with other gauge sectors is anomaly free. Therefore, the anomaly cancellation requirement 
implies that quarks are unified in the following (SU(3)c, SU(3 )l, U(l)x) left- and right-handed representations: 


?L 2 


D 1 ' 2 \ 

( U3 \ 


- U 12 \ : (3, 3*, 0), 

Ql=\D 3 ] : (3, 3,1/3), 

(3) 

. Jl ’ 2 ) L 

\ T Jl 



Dr 2 ’ 3 : (3,1, —1/3), 
4 2 : (3,1,-1/3), 


7yb 2 .3 

U R 

Tr : 


: (3,1,2/3), 
(3,1,2/3). 


(4) 


Here U l L and D l L (i = 1,2,3) are the left handed up- and down-type quarks in the flavor basis. The right handed 
SM quarks U\ l and D l R (i = 1,2,3) and right handed exotic quarks Tr and J R 2 are assigned into SU(3)l singlets 
representations, so that their f7(l)x quantum numbers correspond to their electric charges. 

Furthermore, cancellation of anomalies implies that leptons are grouped in the following (SU(3)c, SU(3 )l, U(1)x) 
left- and right-handed representations: 


r 1,2,3 
h L 


,1,2,3 


e 1 ' 2 ’ 3 

(zH> 2 > 3 'i c , 


e« : ( 1 , 1 ,- 1 ), 

K : ( 1 , 1 , 0 ), 


(1,3,-1/3), 

Pr '■ ( 1 . 1 . — 1 ). 
iV 2 : (1,1,0), 


-ra : (1,1,-1), 

(1,1,0). 


(5) 

( 6 ) 


where v l L and (e/,, tl) are the neutral and charged lepton families, respectively. Let’s note that we assign 
the right-handed leptons as SU(3)l singlets, which implies that their U{l)x quantum numbers correspond to their 
electric charges. The exotic leptons of the model are: three neutral Majorana leptons (t / 1 ’ 2 ’ 3 )^ and three right-handed 
Majorana leptons N R 2 ’ (A recent discussion of double and inverse see-saw neutrino mass generation mechanisms in 
the context of 331 models can be found in Ref. 0 )- 

The scalar sector the 331 models includes: three 3’s irreps of SU(3)l, where one triplet \ gets a TeV scale vaccuum 
expectation value (VEV) v x , that breaks the SU(3)l x U( 1)x symmetry down to SU(2 )l x U(1)y, thus generating 
the masses of non SM fermions and non SM gauge bosons; and two light triplets 7 and p acquiring electroweak scale 
VEVs v v and v p , respectively and thus providing masses for the fermions and gauge bosons of the SM. 

Regarding the scalar sector of the minimal 331 model, we assign the scalar fields in the following [SU(3)l,U(1)x\ 
representations: 


X = 


77 = 


Xi 

X 2 


72 ( Wr 7 4" 4= *Cr7 


V 2 

n° 

V3 


: (3,-1/3), 


: (3,-1/3). 


Pi 


P~ ( ^( v p 4" £p 4= K P ) ) : (3,2/3), 


Pt 


(7) 


We extend the scalar sector of the minimal 331 model by adding the following eleven very heavy SU{3)l scalar 
singlets: 

(1,0), t ~ (1,0), ^ : (1,0), Cj : (1,0), Sj : (1,0), j = 1,2,3. ( 8 ) 

We assign the scalars into T 7 triplet, X 7 antitriplet and T 7 singlet representions. The T 7 (g) Z 2 < 8 > Z 3 <g) Z 14 assignments 
of the scalar fields are: 

7/ ~~ ^lo, 1,^ 3 , 1^ , P ^ ^lo,l, ^ 3 , 1^ , X ~ (1-0,1, 1,1), T — (1 1 , 1, 1, 1) , 

C ~ (3,l,e^,l) , C~ (3,1,1,1), 5-(3,l,e™,l) , a - (l 0 ,l, 1, e"^) . 


(9) 
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It is noteworthy that the SU(3)l singlet scalar field r is the only scalar odd under the Z 2 symmetry and assigned as 
a non trivial T 7 singlet. This scalar field r is crucial for explaining the hierarchy between the SM up and SM down 
type quark masses. 

In the concerning to the lepton sector, we have the following T 7 <g) Z 2 <8> Z 3 <g) Zu assignments: 


Ll ~ ^3,1, e 3 , 1^ , e R ~ ^lo,f, e 3 , , 

( 2ni 2 7 7r \ , „ / 2 tt7 _ \ 

l 2 ,l,e 3 ,e 7 J , jYr ~ (^3,1,e 3 , 1J , 

while the T 7 Z 2 (& Z 3 (& Zu assignments for the quark sector are: 

Ql ~ (lo,l, 1; e 7 ) > Ql ~ 1 ; e 7 ) ) Ql 

^A ~ (V-Re^.e 2 ? 1 ) , (l 0 ,-l,e^,e^) , 

~ (lo,1,1,1), 4 ~ 1) e ^) > 4~ ( 


Hr ~ 3 ,e 7 


~ (lo, 1 , 1 , 1 ), 

4 ~ (lo,l, e 3 ,1 
£>A ~ (lo, — 1 , e 3 
lo,l, 1 , e 7 ) ■ 


)• 

)• 



( 10 ) 


( 11 ) 


Here the dimensions of the T 7 irreducible representations are specified by the numbers in boldface. In the concerning 
to the lepton sector, we recall that the left and right handed leptons are grouped into T 7 triplet and T 7 singlet 
irreducible representations, respectively, whereas the right handed Majorana neutrinos are unified into a X 7 triplet. 
Regarding the quark sector, we assign the quarks fields into trivial T 7 singlet representations. Note that SM right 
handed quarks are the only quark fields transforming non trivially under the Z 3 symmetry. Besides that, let’s note 
that the right handed SM down type quarks are the only fermions odd under the Z 2 symmetry. Furthermore, it 
is worth mentioning that the SU{3)l scalar triplets are assigned to a X 7 trivial singlet representation, whereas the 
SU(3)l scalar singlets are accomodated into two T 7 triplets, one T 7 antitriplet, one T 7 trivial singlet and one T 7 non 
trivial singlet. The SU(3)l scalar singlets X 7 triplets are distinguished by their Z 3 charge assignments. 

With the aforementioned field content of our model, the relevant quark and lepton Yukawa terms invariant under 
the group Q, take the form: 


AQ) 

l-y 


vVPQ^ukj 4 + vWq'lp'Urj 8 + yii ] Q 2 Lp* ul R^3 + v™QIp*Ur^2 


+v% ) Q 1 lp*u 3 R j 2 + viVQlyU 1 ^ + + y^QlnUl'- 


X u )n 3 


(U)-. 


Z. 


+ 2 / 334 l ? ?4 + y (1 ^lXTr + y[ 'QlX*Jr+V2 ’QlX*Jr + 3/33 'QlPDr 


■ XT) 


,,(/) 


* t1 


■X J )1 


* t2 


,(-0)7 


A 3 


+y[? ) Q 1 LV*D 1 R ^ + v^qWdI^- + 


3_3 


A 7 


+Z/13 ^QlV*D r —~s —b 2/37 ’Q L pD R —-g —b y 23 ’Ql^Dr—t —b 2/32 ’QlPDr— 7—b H.c, 


A 6 

_2_3 


A 6 


+ y i 22 ) Q 2 L v*D 2 R - 


A 5 


(0)7 


CFT 


A 5 


A 5 


A 


A 


( 12 ) 


~C[ L) = h%> (L L p£) lo e R ji + h$ (L L tf) u ii R ji + h$ (L lP ()^t r ^ 

+h^ ( L l xN r ) ±o + -hi N ( N r N r ) 3 £ + h 2N (AT r N r )^S 

+h p e abc (L a L (L c L ) b ) 3 ( P *) c j+H.c, (13) 

where y^’ D ^ ( i,j = 1,2,3), h p e\ h p JA, h P r\ h^\ h\N , h 2 N and h p are 0(1) dimensionless couplings. 

In the following we explain the role each discrete group factors of our model. The X 7 and Z 3 discrete groups 
reduce the number of the SU(3)c < 8 > SU(3)l < 8 > U(l)x model parameters. This allow us to get predictive and viable 
textures for the fermion sector that successfully describe the prevailing pattern of fermion masses and mixings, as 
we will show in sections EH and m We use T 7 since it is the minimal non-Abelian discrete group having a complex 
triplet 0 , where the three fermion generations can be naturally unified. The Z 3 symmetry determines the allowed 
entries of the neutrino mass matrix. Furthermore, the Z 3 symmetry distinguishes the right handed exotic quaks, 
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being neutral under from the right handed SM quarks, charged under this symmetry. This results in the absence 
of mixing between SM quarks and exotic quarks. Consequently, the symmetry is crucial for decoupling the SM 
quarks from the exotic quarks. The Z -2 symmetry separates the right handed SM down type quarks, odd under 
this symmetry, from the SM up type quarks, even under Z 2 . Consequently, the Z 2 symmetry is responsible for the 
mass hierarchy between SM up and SM down type quarks. Note that the SU(3)l scalar singlet r, the only X 7 
non trivial singlet and Z 2 odd scalar, only appears in the SM down type quark Yukawa terms. The Z14 symmetry 
generates the hierarchy among charged fermion masses and quark mixing angles that yields the observed charged 
fermion mass and quark mixing pattern. Furthermore, it is noteworthy that the five dimensional Yukawa operators 
X (TlpC) i 0 e Ri x (TlpC) i 2 P r an d x; {Llp£) 1i t R are invariant under T 7 but not under the Z14 symmetry, because 
the right handed charged leptons are Z\± charged. We use Z\± because it is the smallest lowest cyclic symmetry, 
that allows to build a twelve dimensional charged lepton Yukawa term, from a ^7 insertion on the ^ eu 

operator. That aforementioned twelve dimensional charged lepton Yukawa term is crucial to explain the smallness of 
the electron mass, without tuning its corresponding Yukawa coupling. 

To get a predictive model that successfully accounts for fermion masses and mixings, we assume that the SU (3 )l 
singlet scalars have the following VEV pattern: 

(a)=v a e i *, (t) = v T , <0 = ^| (1,1,1), (C) = ^| ( 1 , 0 , 1 ), (S) = ^ |(1,1,-1). (14) 

We justify this choice of directions in the T 7 space by the observation that they describe a natural solution of the scalar 
potential minimization equations. Indeed, in the single-field case, T 7 invariance readily favors the (1,1,1) direction 
over e.g. the (1,0,0) solution for large regions of parameter space. The vacuum (£) is a configuration that preserves 
a Z 3 subgroup of T 7 , which has been extensively studied by many authors, (see for example Ref [35lp. In the next 
subsection, we consider the minimization conditions of the high energy scalar potential C3 of our model, and show 
that our chosen VEV directions for the two X 7 triplets, i.e., £, S and the T 7 antitriplet £ scalars in Eq. m , are 
consistent with a global minimum of this scalar potential. 

Besides that, the SU(3)l scalar singlets are assumed to acquire vacuum expectation values at a very high energy 
Am* ^ v x ss 0(1) TeV, excepting £ • (j = 1,2,3), whose vacuum expectation value is much lower than the scale of 
electroweak symmetry breaking v = 246 GeV. Let’s note that at the scale A int , the SU (3)c < 8 > SU (3)^ ® U (l)x ®Ty® 
Z 2 <8> Z 3 (§) Z 14 symmetry is broken to SU(3)c <8> SU(3)l <S> U( 1)x by the vacuum expectation values of the SU(3)l 
singlet scalar fields £ •, Sj, cr and r. 

Considering that the charged fermion mass and quark mixing pattern arises from the Z 2 < 8 > Z 3 ® Z \4 symmetry 
breaking, we set the VEVs of the SU(3)l singlet scalars S, £, a and r, as follows: 

VS ~ V£ = v a = V T = A int = AA, (15) 

being A = 0.225 one of the Wolfenstein parameters and A our model cutoff. Consequently, the VEVs of the scalars in 
our model have the following hierarchy: 


V( « v p ~ v v ~ v « v x « Am t . (16) 

Thus, the SU(3)l scalar singlets having Yukawa interactions with the right handed Majorana neutrinos get VEVs at 
very high scale, then providing very large masses to these Majorana neutrinos, and thus giving rise to a double seesaw 
mechanism of active neutrino masses. Consequently, the neutrino spectrum includes very light active neutrinos as 
well as heavy and very heavy sterile neutrinos. As we will shown in detail in the next section, the smallness of the 
active neutrino masses is attributed to their scaling with inverse powers of the high energy cutoff A as well as by their 
quadratic dependence on the very small VEV of the Z 2 <Z> Z 3 ® Zu neutral, SU{3)l singlet and T 7 antitriplet scalar 
field (■ 


B. High energy scalar potential 


As previously mentioned, all singlet scalars, excepting £ • (j = 1, 2, 3), acquire vaccum expectation values, much larger 
than v x , which implies that the singlet scalars are very heavy and thus the mixing between these scalar singlets and 
the SU ( 3) l scalar triplets can be neglected as done in Ref. 3C|. For simplicity we assume a CP invariant scalar 
potential with only real couplings as done in Refs. [HI, [ 2 ^, [30j, [42l- The high energy scalar potential which involves 
only the SU ( 3) L singlet slcalars is given by: 
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Vi = 4(ti*) l0 +»s(SS*) l0 +»H CC*) 1o +/4DD + /4(tt*) 

+ pi (cd 3 c+ c 2 (ss*) 3 c+c 3 (^) 3 c+c 4 (co 3 r+c 5 (so 3 s*+c 6 (r o 3 c 

+ c 7 (S*0 3 S + c 8 (ec)a S + Cg (SC) 3 C + Cio (CC) 3 C + D (CO3C* + Do (CC*k C 

+ c 16 (cn 3 c*+c 12 (ts) 3 c +C13 (ro 3 s*+c 14 (s* o 3 c* + +«i (cd 1o (cc*)i 0 

+« 2 (enu (cc*) l2 +« 3 (cd 3 (cd 3 +m (55*) lo (ss*) l0 + « 5 (^*) 11 (ss*) l2 
+« 6 (55*) 3 (55^3+« 7 (cn l0 (cn l0 +« 8 (cc*) 4l (co l2 +« 9 (cn 3 (cn 3 
+[« 10 (co 3 (cck+775 (6 s) 3 (Dk+7i (co 3 (ss ) 3 + 72 (ss) 3 (cck + ha 

+73 (CC*) lo (S5*)i 0 + 74 (CC*)^ (SS*) l2 + 75 (CD la (SS*) U + 76 (CDs (SS*) 3 

+77 (cc*k (ss*) 3 + [714 (co 3 (ss) 3+ 7l5 (co 3 (cck+ ha + 7 8 (d*)i 0 (sd 1o 

+79 (m u (sn l2 +7h (d*) 12 (5n Xl +7 i 2 m a (cs*) 3 + 7 i 3 (c*d 3 (c% 

+717 (D*) 3 (CS) 3 + 718 (C* 5)3 (£S*k + 719 (CDi 0 (COi 0 + 720 (CDi t (CC*)i 3 
+721 (CDi a (CC*)i x + 723 (CDs (CC*) 3 + 724 (CDs (CC*k 
+ [725 (CC *) 3 (CO3 + 726 (CC ) 3 (CC *) 3 + HA + 730 (CC ) 3 (C*C*k 

+ 72 r (co 10 (c*d 1o +728 (ccd (c*c*) l2 + 7 29 (co l2 (ctd 
+ 731 (tn 3 (cck+[732 (cck (c*ck+733 (c*o 3 (cck + ha 

+734 (SS*) 1 0 (CC*)l 0 + 735 (SS*) U (CD 1 2 + 736 (SS*) 1 , (CDu + 737 (SS*) 3 (CC*k 
+738 (CC*k (S5 *) 5 + [739 (SS*) 3 (CCk + 740 (CC ) 3 (55*)a + JT.c] + 744 (Dk (S*C*k 
+741 (D)l 0 (S*C*) lo + 742 (SQ ll (S*C)i 2 + 743 ( s 0i a (S*C) U + 745 D*C*k (Dk 
+ [749 (Dk (CC*k + 750 (CDs (D) 3 + 751 (D ) 3 (CCk + 752 (CC ) 8 (D)s + JY.c] 

+ [753 (CC) lo (SQlo + 754 (CC)u (SC) l2 + 755 (CC) la (SOu 
+ 756 (CC) 3 (SO 3 + 757 (D) 3 (CCk + -ff-c] + 766 (CDl 0 DD + 7l6 (CC*)l 0 (°°*) 

+758 (CDl 0 (™*) + 759 (S5*) l0 (o-o-*) + 716 (CDl, C 2 + 767 (CC*)l 2 D 

+761 (CC*)l 0 (TT*) + 762 (SS*) 1 0 (TT*) + 7 68 (crcr*) 2 + 769 ( TT *f + 770 DD ( TT *) 

+771 (CDi, r 2 + 772 (CC*)i 2 r* 2 + 773 (55*) Xl r 2 + 774 (SS*) l2 r* 2 (17) 

Now we are going to determine the conditions under which the VEV pattern for the SU(3)l singlet scalars, given in 
Eq. (1141) . is a solution of the high energy scalar potential. In view of the very large number of parameters of the high 
energy scalar potential, in order to simplify the analysis, we assume universality in its trilinear and quartic couplings, 

i.e. 


Ki = 7 j = n, i = (1 - 10), j = (1 - 75), 

C k = C, k = (1 -15) (18) 


Then, from the minimization conditions of the scalar potential and taking into account our assumption that v^ j 
(j = 1, 2, 3) are much lower than the electroweak symmetry breaking scale v = 246 GeV, the following relations are 
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obtained: 

d{Vi) 

dv 


d(Vi) 

dv t 2 


d(Vi) 

dv H 3 


8 Kvf ± + Ahv\ 2 + 2v^/il + 2kv £l (yl 2 + vfj + 12 «u| i ^ 3 + 2 /m^v^ + 2/cv^ (7v Sl + 3u ?3 ) 

+2ct| 2 (2n ?1 + 3n ?3 ) + 4/cvf^ + 2 k + 3^) u 5l + 4 ku Si (n ?2 t;s 2 + ^ 3 «s 3 + 3u €l us 3 ) 

0 , 

8 K y l a + 4 kv^ + 2v£ 2 fj, 2 + 2 K v i2 (v 2 ^ + vf 3 ) + 12 kv^v (i + 2kv^ + 2ku| 2 (7v ? 2 + 3t> Cl ) 

+ 2 kv | 3 (2n ?2 + 3n ?3 ) + 4 kv| 1 ^ 2 + 2ki^ 2 + 4kus 2 (ug^Sj. + v ?3 us 3 + 3u| 2 w Sl ) 

0 , 

8 Kn | 3 + 4 kv| i + 2 v ^ n \ + 2 k ^ 3 + v | 2 ) + 12 k ^| 3 ^ 2 + 2 k ^ 3 «| 2 + 2 kv § 3 ( 7 v ?3 + 3 u ?2 ) 

+2 kw|i (2 n ?3 + 3n ?1 ) + 4 ku| 2 ^ 3 + 2 k^v ?3 + 4kus 3 (u$ 2 t;s 3 + + 3W£ 3 us 2 ) 

0 , 


(19) 


W> 

9^Si 


W) 

9vs 2 


d(Vi) 

dvs 3 


8kv 3 Si + 4kw| 2 + 2v Sl n 2 + 2 kv Si (vg 2 + Vg 3 ) + 12 kv 2 Si vs 3 + 2nv Sl Vg 2 + 2 kv^ ( 7v Sl + 3us 3 ) 

+2ct| 2 (2v Sl + 3ws 3 ) + 4ku| 3 U Si + 2k (v 2 a + 3v 2 ) v Sl + 4ku 5i ( vg 2 v i2 + vg 3 v i;i + Sv Sl v i3 ) 

0, 

8ki>s 2 + 4kv^ 3 + 2 vs 2 n 2 + 2kvs 2 + i'g 3 ) + 12kv^ 2 vs 1 + 2nvs 2 Vg 3 + 2kv^ 2 (7 vs 2 + 31^) 

+2 kv 2 3 (2 v S2 + 3 v S3 ) + 4kv\_v S2 + 2nv 2 vg 2 + 4kd $2 (v Sl v il + vg 3 v i3 + 3v S2 v^) 

0, 

8kv 3 S3 + 4kv 3 Si + 2vs 3 /4 + 2kvs 3 (v 2 Sl + u| 2 ) + 12nv 2 g 3 v S2 + 2nvs 3 v 2 g 2 + 2kv\ 3 (7 v Ss + 3vg 2 ) 

+2ku| i (2 vs 3 + 3v Sl ) + 4Kw| 2 n5 3 + 2 kv 2 vs 3 + 4 kv^ 3 (vs 2 v^ 2 + vs^ + 3 vs 3 v^ 2 ) 

0, (20) 


d(Vi) 

dv (l 

d(Vi) 

dv <2 

d(Vi) 

dv C 3 


6 C + v Sl v S2 + 2 (v£ 3 vs 2 + v^ 2 vs 3 ) + ^ 3 vs 3 ] + 2 ^ 2 v Cl = 0 , 

6 C [v| 2 Uf 3 + vs 2 vs 3 + 2 (v£ 3 v Sl + V^VSa) +v^v Sl ] + 2 /j, 2 v ^ 2 = 0 , 

6 C + v Sl vs 3 + 2 (v£ 2 v Sl + v^v S2 ) +Wf 2 t>s 2 ] + 2 Mc v c 3 = °> 


( 21 ) 


d(Vi) 

dvfj 

d(Vi) 

dv T 


2 v a 

2v t 


+ K (v 2 ^ + w | 2 + vl 3 + v 2 Si + v 2 s2 + v 2 g 3 + 2 v 2 a + v 2 ^ 


= 0 , 


( 3v l 


+ 3v Sl +v a + 2v 2 


= 0 . 


( 22 ) 


From the expressions given above, and using the vacuum configuration for the SU(3)l singlet scalars given in Eq. 
m , we find the following relations: 

(60u| + 41r)|) , ns = ~^ K ( 59v ( + 10 v S ) , 

/*C = _—^-( v e+ v s), nl = Mr = ~k(4v\ +4) > 

2 v 2 — 2 Vg + v^vg = 0 . 


( 23 ) 


















Taking the positive solution of the previous equation, we find: 


v s 


1 + ^5 

— 


~ 0.81u^, 


(24) 


which is consistent with our previous assumption described by Eq. (USD- Our results show that the VEV directions 
for the two T 7 triplets, i.e., £, S and the T 7 antitriplet £ scalars in Eq. (ED. are consistent with a global minimum of 
the high scalar potential ED of our model, for a not fine-tuned region of parameter space. 


C. Low energy scalar potential 


The renormalizable low energy scalar potential of the model takes the form: 

V H = /4(X f X) + + Mp(pV) + / (: XiVjPk^ 3k + H.c) + Ai(x t x)(x t x) 

+A2 (pVXpV) + A 3 (?? t ? 7 )(r? t ? 7 ) + A 4 (x t x)(p t p) + ^(x'xKv^v) 

+A 6 (p t p)(r 7 t 77 ) + A 7 (x t ? 7 )(? 7 t x) + A 8 (x t />)(P t x) + Mp^v)(v*p)- (25) 


After the symmetry breaking, it is found that the scalar mass eigenstates are connected with the weak scalar states 
by the following relations: [l4 . 0 : 


with 



(26) 

(27) 


_( cosa T ((3 T ) sma T (P T )\ n -(~ 1 

“t(/ 3 t ) ^ — gin a T ( P T ) cosq;t(/3t) / 0 1 J ’ 

where tan/3 T = v^/vp, and tan 2 olt = Mi/(M 2 — M 3 ) with: 

Mi = 4A QV^Vp + 2 V 2 fv x , 

Mi = 4:\ 2 Vp - V2fv x t&n/3 T , 

M 3 = 4A 3 v^ — V 2 fv x /tan/3 T . 


(28) 


(29) 


The low energy physical scalar spectrum of our model includes: 4 massive charged Higgs (ff±, H 2 ), one CP-odd Higgs 
(Aq), 3 neutral CP-even Higgs (h°, H®, H®) and 2 neutral Higgs (H 2 , H 2 ) bosons. The scalar h° is identified with the 
SM-like 126 GeV Higgs boson found at the LHC. It it noteworthy that the neutral Goldstone bosons G?, G(j, G 2 , 
G 2 are associated to the longitudinal components of the Z, Z', K° and K gauge bosons, respectively. Furthermore, 
the charged Goldstone bosons Gf and Gf are associated to the longitudinal components of the W ± and K ± gauge 
bosons, respectively jlH fTbf . 


III. LEPTON MASSES AND MIXINGS 


From Eqs. ED. Ell . (El) and using the product rules of the T 7 group given in Appendix [A] it follows that the mass 
matrix for charged leptons is: 


Mi = V^Pidiag (m e , m M , m T ) 


VlL y/3 


/I 

1 

1 \ 

, j 

/ e 7i0 

0 

0 \ 

1 

UJ 

w2 

’ Pl = 

0 

e a<p 

0 

VI 

CO 2 

- 

1 1 

l 0 

0 

e 2 ^ / 


u> = e 3 


(30) 
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where the charged lepton masses read: 


m e = h$\ s ^ 


= h^X 5 ^- 
«* v^’ 


m T - hfiX 3 ^. 


(31) 


Taking into account that v p ~ v = 246 GeV, it follows that the charged lepton masses are related with the electroweak 
symmetry breaking scale by their scalings with powers of the Wolfenstein parameter A = 0.225, with 0(1) coefficients. 
In the concerning to the neutrino sector, the following neutrino mass terms arise: 



(32) 


where the TV discrete flavor group constrains the neutrino mass matrix to be of the form: 

hpVpV f 



CO 

X 

CO 

0 

Mr, 

03x3 

\ 

M v = 

II 

03x3 

M x 



\ 03x3 

K 

Mr 

) 


h V * 

~ ftl»-7= 

{ 1 

— X 

X 

Mr - 

-X 

V X 

1 

X 

X 

1 


M d = 



x = 


h-2NVS 
hl N V£ ' 


M * = h * L) f2 



(33) 


Since the SU(3)l singlet scalars having Yukawa interactions with the right handed Majorana neutrinos acquire VEVs 
at very high scale, these Majorana neutrinos are very heavy, so that the active neutrinos get small masses via a double 
seesaw mechanism. 

The full rotation matrix, which diagonalizes the neutrino mass matrix, takes the approximate form p~ 0 ] : 


u = 

/ 14 b 2 u x 0 \ 

-B\v v U x BiUr , 

(34) 


V 0 B\u x Ur ) 


where 

h ta_ 

II 

£ 

& 1 

*4 

b\ = m d {m%) 1 MrM~\ 

(35) 

and the neutrino mass matrices for the physical states are: 

MW = Md(mT)' 1 MrM-'mZ, 

(36) 

Mf> = 

-M x M^M^, 

(37) 

m (3) = 

Mr , 

(38) 


being the mass matrix for light active neutrinos, while MI) 2 ' 1 and M^ are the heavy and very heavy sterile 

neutrino mass matrices, respectively. Consequently, the double seesaw mechanism gives rise to light active neutrinos 
as well as to heavy and very heavy sterile neutrinos. Moreover, the neutrino mass matrices M^ and Ml 3 ' 1 are 

diagonalized by the rotation matrices 14, Ur and U x , respectively. 141]. 

Using Eq. (1361) . we find the following mass matrix for light active neutrinos: 


where 


(A 0 A\ 
Ml l) = 050 , 
\A 0 A J 


h 1 N h 2 p v%v%V£ 
2sJ?>h { x ' ) v\ A 2 ’ 


B = 




(JliNV£ + Jl2NVs) ■ 


(39) 


(40) 


From Eq. (13*71) it follows that the light active neutrino mass matrix only depends on two effective parameters: A and 
5, which determine the neutrino mass squared splittings. Let’s note that A and B are supressed by their scaling with 
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inverse powers of the high energy cutoff A. Furthermore, we have that the smallness of the active neutrino masses 
arises from their scaling with inverse powers of the high energy cutoff A as well as from their quadratic dependence 
on the very small VEV of the Z 2 ®Z%® Z 14 neutral, SU (3 )l singlet and T 7 antitriplet scalar field Considering that 
the orders of magnitude of the SM particles and new physics yield the constraints v x > 1 TeV and u 2 + v 2 p = v 2 and 
taking into account our assumption that the dimensionless lepton Yukawa couplings are 0(1) parameters, from Eq. 
(l39l) and the relations = AA, v p ~ 100 GeV, v x ~ 1 TeV, we get that the mass scale for the light active neutrinos 

v 2 

satisfies m v ~ 10~ 3 -^-. Consequently, setting = 1 GeV, we find for the cutoff of our model the estimate 

A ~ 10 5 TeV, (41) 

which is of the same order of magnitude of the cutoff of our S 3 lepton flavor 331 model |2{J. Consequenty, we find that 
the heavy and very heavy sterile neutrinos have masses at the ~ MeV and ~ TeV scales, respectively. Furthermore, 
from the aforementioned considerations, as well as from Eqs (fl5l) . (f35ll and ED. it follows that: 




v p v c v <; ~ 10 —5 


<4 A 


<< 


(Bi)i 


~ 10 " 3 , 

V Z 


i,j = 1,2,3, 


(42) 


Moreover, we hnd that the mass matrix for light active neutrinos is diagonalized by a rotation matrix V„, as 
follows: 


V V 1V± V V v 


( mi 

0 

0 \ 


( cos 9 

0 

sin 9 \ 

° 

m 2 

° ’ 

with V v = 

° 

1 

° 

V 0 

0 

m 3 / 


V — sin 9 

0 

COS0 / 


-4 


(43) 


where 0 = + 7 r /4 and 6 = — 7 r /4 correspond to normal (NH) and inverted (IH) mass hierarchies, respectively. The 
masses for the light active neutrinos, in the cases of normal (NH) and inverted (IH) mass hierarchies, read: 


NH : 

: 0 = 

7 r 

= + 4 : 

m Vl = 0, 

m„ 2 = B, 

TYii /3 — 2^4., 

(44) 

IH : 

: 9 = 

7 r 

~4 ’ 

TTLi/^ — 2A, 

m „ 2 = B , 

m U3 = 0. 

(45) 


Besides that, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix has the following form: 


u = V? L P t V v = 


e 7i * cos 8 

e 2t ^ sin 6 

e 4 ^ 

2 zd> rt 

e ^ cos 6 1 

e 7tt ^ sin 6 

V3 

V3 

Vs 

V3 


V3 

e 7x< ^ cos 6 

e 2 ^ +2 r sin g 


e 2 ^+f 1 

cos 6 

, e 7i * sin g 

Vs 

Vs 

Vs 

Vs 


1 Vs 

e 7x< ^ cos 6 

e^-^sin g 

e 4 ^ +2 r 


cos 0 

1 e 7l< ^ sin 0 

Vs 

Vs 

Vs 

V3 


1 Vs 


(46) 


Note that while the PMNS leptonic mixing matrix only depends on a single parameter 0 (which arises from the 
charged lepton mass matrix), the neutrino mass squared splittings are determined by two parameters, i.e., A and B. 
Furthermore, we find that the lepton mixing angles are given by: 


sin 2 6*12 


\Ue 2 \ 2 = 1 

1 — |C/es| 2 2 T cos 50’ 


sin 2 9 13 
sin 2 023 


\Ue3\ 2 = ^(1 ± COS 50), 

1 Ufj. 3 1 2 _ 1 ± V3 sin 50) 

1 — | U e 3 1 ^ 2 4 q= 2 cos 50' 


(47) 

(48) 

(49) 


Then, from Eq. it follows that the limit 0 = 0 and 0 = 7r for the inverted and normal neutrino mass hierarchies, 
respectively, correspond to the trimaximal mixing, which predicts a vanishing reactor mixing angle. Let’s note that 
the mixing angles for the lepton sector only depend on a single parameter (0), while the neutrino mass squared 
splittings are controlled by two parameters, i.e., A and B. 
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The Jarlskog invariant and the CP violating phase are given by 


J = Im {u el u,aU^u; A ) = 




sin 5 = 


_8J_ 

cos 0 i 3 sin 2012 sin 20 2 3 sin 2013 ' 


(50) 


From the relation 0 = ±^, we predict J = 0 and 5 = 0, which implies a vanishing leptonic Dirac CP violating phase. 

In the following the three free effective parameters 0, A and B of the SM lepton sector of our model are adjusted to 
accommodate the experimental values of three leptonic mixing parameters and two neutrino mass squared splittings, 
shown in Tables mini for the normal (NH) and inverted (IH) neutrino mass hierarchies, respectively. The parameter 
<f> is adjusted to reproduce the experimental values of the leptonic mixing parameters sin 2 6 ij , whereas A and B for 
the normal (NH) and inverted (IH) neutrino mass hierarchies read: 


NH : 

m Vl = 0, 

m„ 2 = B = yj Am^ « 9meV, 

3 

w 

II 

to 

£ 

II 

J Am 31 S 

s 50meV; 

(51) 

IH : 

: m„ 2 = B 

= J Am 2 ! + A m 2 3 k, 50meV, 

II 

to 

A 

II 

J Ato 13 ? 

a 49meV, m V3 = 0, 

(52) 


which follows from Eqs. (l45l) . (l44l) and the definition A= mf — m 2 . We take the best fit values of Am^ from 
Tables H and El for the normal and inverted neutrino mass hierarchies, respectively. 

To reproduce the experimental values of the leptonic mixing parameters sin 2 0 ly given in Tables id hd we vary the 
<j> parameter, finding the following result: 

NH : cf) = 0.576 7r, sin 2 0 i2 « 0.34, sin 2 0 23 « 0.61, sin 2 0 i3 « 0.0232; (53) 

IH : (t>= 0.376tt, sin 2 0 i2 « 0.34, sin 2 0 23 « 0.61, sin 2 0 i3 « 0.0238. (54) 


Consequently, we find that sin 2 0i3 is in excellent agreement with the experimental data, for both normal and 
inverted neutrino mass hierarchies, whereas sin 2 0i 2 and sin 2 0 2 3 stay in the experimentally allowed 2 a range. Thus, 
our predictions for the neutrino mass squared splittings and leptonic mixing parameters, are in very good agreement 
with the experimental data on neutrino oscillations, for both normal and inverted mass hierarchies. Furthermore, 
another relevant prediction of our model is a vanishing leptonic Dirac CP violating phase. 


Parameter 

Am| 1 (10“ 5 eV 2 ) 

Am| 1 (10- 3 eV 2 ) 

(sin 2 0 12 ) exD 

(sin 2 0 23 ) exD 

(sin 2 0 i 3 ) bxd 

Best fit 

7.60 

2.48 

0.323 

0.567 

0.0234 

lcr range 

7.42 - 7.79 

2.41 - 2.53 

0.307 - 0.339 

0.439 - 0.599 

0.0214 - 0.0254 

2 cr range 

7.26 - 7.99 

2.35 - 2.59 

0.292 - 0.357 

0.413 - 0.623 

0.0195 - 0.0274 

3cr range 

7.11 - 8.11 

2.30 - 2.65 

0.278 - 0.375 

0.392 - 0.643 

0.0183 - 0.0297 


Table I: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from Ref. 
[25| . for the case of normal hierarchy. 


Parameter 

Am 21 (10 5 eV 2 ) 

Amf 3 (10 3 eV 2 ) 

(sin 2 0i 2 ) 

V / exp 

(sin 2 0 23 ) exD 

( sin2 0 13)exo 

Best fit 

7.60 

2.38 

0.323 

0.573 

0.0240 

lcr range 

7.42 - 7.79 

2.32 - 2.43 

0.307 - 0.339 

0.530 - 0.598 

0.0221 - 0.0259 

2 cr range 

7.26 - 7.99 

2.26 - 2.48 

0.292 - 0.357 

0.432 - 0.621 

0.0202 - 0.0278 

3cr range 

7.11 - 8.11 

2.20 - 2.54 

0.278 - 0.375 

0.403 - 0.640 

0.0183 - 0.0297 


Table II: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from Ref. 
[25j, for the case of inverted hierarchy. 


Now we determine the effective Majorana neutrino mass parameter, which is proportional to the neutrinoless double 
beta (Ov/3/3) decay amplitude. The effective Majorana neutrino mass parameter is given by: 

y ' U e k m v k i 
j 


rripp = 


(55) 
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where [/A and m^ k are the PMNS mixing matrix elements and the Majorana neutrino masses, respectively. 

From Eqs. flU, m , (El and El- we find that that the effective Majorana neutrino mass parameter, for both 
Normal and Inverted hierarchies, takes the following values: 


TO/3/3 


3 meV for NH 
40 meV for IH 


(56) 


Our obtained values mpp s=s 3 meV and rnpp ss 40 meV for the effective Majorana neutrino mass parameter, for normal 
and inverted hierarchies, respectively, are beyond the reach of the present and forthcoming Ou/3/3 decay experiments. 
The current best upper bound on the effective neutrino mass is mpp < 160 meV, which corresponds to Ty^ ( 136 Xe) > 
1.6 x 10 25 yr at 90% C.L, as indicated by the EXO-200 experiment [iij . This bound will be improved within a not 
too far future. The GERDA “phase-II” experiment [45], [46[ is expected to reach Ty^( 7e Ge) > 2 x 10 26 yr, which 
corresponds to mpp < 100 meV. A bolometric CUORE experiment, using 130 Te [13|, is currently under construction 
and has an estimated sensitivity of about T°^ 2 /3/3 ( 130 Te) ~ 10 26 yr, which corresponds to 171 / 3/3 < 50 meV. Furthermore, 
there are proposals for ton-scale next-to-next generation 0 v/3f3 experiments with 136 Xe [3, [49] and 76 Ge mutt) 
claiming sensitivities over Ty / 33 ~ 10 27 yr, which corresponds to mpp ~ 12 — 30 meV. For a recent review, see for 
example Ref. [ 5 ]]]. Consequently, as follows from Eq. (l5f>ll . our model predicts Ty 33 at the level of sensitivities of 
the next generation or next-to-next generation 0^/3/? experiments. 

Regarding the sterile neutrino sector, from Eqs. El and (13811 we find that the sterile neutrino mass matrices are 
given by: 


M, (2) = - 


= hlN ^L 



X = 


: — 1 


2 x 2 + x — 1 : 


Y = 


2 x 2 + x — 1 


(57) 


(58) 


(2) (3) 

The sterile neutrino mass matrices A11 and Ml are diagonalized by a rotation matrix Ur = U x , according to: 

k = 2,3, (59) 



( m[ k) 

0 

0 \ 


/ 

1 

y/s 

1 

V2 

V6 ^ 

UrM^Ur = 

0 

m 2 (k) 

0 

II 

II 


1 

y/3 

1 

V2 

1 

VE 


l 0 

0 

a 4 k) ) 


l 

1 

Vs 

0 

-,/I 

V 3 / 


where the sterile neutrino masses are given by: 

. 2 


M 1 (2) = 


m} 3) = 


^(hpy 


v: 


2hi N V£ 

hi (1 - 2 x) 


■ (X - 2 Y) 


M. 


(2) 


= iff' = 


V3(h[ L) )‘ 


2h\NV£ 


' {X + Y) 


m 2 (3) = m 3 (3) = 


hiN— 7 = (1 + x) 


Vz 


(60) 

(61) 


which implies that the heavy and very heavy sterile neutrino spectrum includes two degenerates heavy neutrino and 
one light neutrino states. 


Furthermore, as follows from Eq. El and the relation (_B 2 )iJ ~ 10 5 << (F?i)^ ~ 10 3 (i,j = 1,2,3) given by 

Eq. HUD, we can connect the neutrino fields vl = {v\ l, ^ 2 l, l) T , v r = (vi R , v ( ,y R . v^r) and Nr = (N^ r , N^r, N^ r ) 
with the neutrino mass eigenstates by the following approximate relations: 



( 

,,c 


V R 

- 

\ N K) 

V 




Ur d 2) + BiUrZ {3) 


>L 

t(3) 


L 

( 2 ) 


Aj) _ 

SL — 


/ Aj) 

51 L 

AA 

52 L 


3 = 1,2,3. 


(62) 





























13 


where 47 4i and 4 l (k = 1, 2, 3) are the light active, heavy sterile and very heavy sterile neutrinos, respectively. 
As previously mentioned, the heavy sterile neutrinos have MeV scale masses and thus correspond to dark matter 

(3) 

candidates. Furthermore, we assume that the lightest of the very heavy sterile neutrinos, i.e., Tj has a TeV scale 
mass and thus corresponds to a candidate for detection at the LHC. 


IV. QUARK MASSES AND MIXING. 


From Eq. m and taking into account that the VEV pattern of the SU (3 ) L singlet scalar fields is described by Eq. 
m , with the nonvanishing VEVs set to be equal to AA (being A the cutoff of our model) as shown in Eq. (1151) . it 
follows that the SM quark mass matrices have the form: 



/ 

47 a 4 

47a 3 

47 a 2 ^ 

/o’ 



47 a 7 

47 a 6 

47a 5 

M v = 


47 a 3 

47 a 2 

47a 

m d = 


47a 6 

47a 5 

47 a 4 


V 

47a 2 

47a 

47 J 

V 2 


\ 

47a 5 

47a 4 

47 a 3 


v 

71 ’ 


(63) 


where A = 0.225 is one of the Wolfenstein parameters, v = 246 GeV the scale of electroweak symmetry breaking and 
(*i j = 1,2,3) are 0(1) parameters given by the following relations: 


a. 


(U) 

nj 


— V P J(G-i-n) 

i/n? c 
J v 


n (°) _ 

- y nj y e ’ 


rij 


JU) _ i(3-j)4> 

Z 3 j ~ y e > 

S D ) _ 

a 3j ~ 2/3 j „ e 


n = 1,2, 


3 = 1,2,3. 


(64) 


Furthermore, the exotic quark masses read: 

,.(T) V X 


m T = y ' 


V2’ 


( J ) v x 
mji=y , 


(/) 

y\ 


y 


(T) 


m T , 


m — -A/) Vx 


( J ) 
2/2 


y 


(T) 


rriT- 


(65) 


Since the charged fermion mass and quark mixing pattern arises from the breaking of the Z 2 ® ® Z\± discrete 

group, we asssume an approximate universality in the dimensionless SM quark Yukawa couplings, as follows: 


W = a m p i 0i 


SU) 


= a. 


( u ) 


x (C/) 

x 33 


= a 


(U) 


i (U) 


z (D) 

ni 


= a 


(U) 


% {D) 

l 2'2 


1- 


i-l! 

2 


SU) 

i 13 


= a 


(U) 


1- 


a {D) ~ a (D) 

u 23 — u 33 


1 -- 

2 


(S u ) _ 

a 23 — 


(U, D ) (U,D) 


1 {U) 

‘13 


1 - 


i,j = 1,2,3, (66) 


It is noteworthy that exact universality in the dimensionless quark Yukawa couplings predicts massless up, down, 
strange and charm quarks, in contradiction with the experimental data on quark masses. Consequently, to generate 
these masses, the universality in the quark Yukawa couplings has to be broken. Besides that, for simplicity, we assume 
that the complex phase responsible for CP violation in the quark sector only arises from up type quark Yukawa terms, 
as indicated by Eq. (1661) . In addition, for the sake of simplicity, we fix a ’*/' 1 = 1, which is suggested by naturalness 
arguments. Therefore, the SM quark mass matrices are given by: 



( 68 ) 
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Let’s recall that the quark sector has 10 effective parameters, i.e, A, a^\ a( T \ a^, a^\ a [^, and the 

phases fa and fa to describe the quark mass and mixing pattern, which is determined by 10 observables. Nevertheless, 
not all these effective parameters are free since the parameters A and o^' 1 are fixed while the remaining 8 parameters 
are adjusted to reproduce the physical observables in the quark sector, i.e., 6 quark masses and 4 quark mixing 
parameters. The results shown in Table Hill correspond to the following best-fit values: 

a{ J) ~ 0.64, af ] ~ 0.77, fa ~ -9.03°, fa ~ -4.53°, 

a[ D 2 ] ~ 2.03, a ~ 1.75, a[^ ~ 1.15, affi ~ 1.40. (69) 


Observable 

Model value 

Experimental value 

m u (MeV) 

1.59 

1 AK+0.56 

0.45 

rric(MeV) 

673 

635 ± 86 

m t (GeV) 

180 

172.1 ±0.6 ±0.9 

nid(MeV) 

2.9 

2 Q + 0 - 5 
z<y -0.4 

m a {MeV) 

59.7 

C 7 7 +I 6.8 
°* *'-15.7 

m b {GeV ) 

2.98 

9 29+0.09 

Z.OZ—0,04 

|Wd| 

0.975 

0.97427 ± 0.00015 

Ks 

0.224 

0.22534 ±0.00065 

\v ub \ 

0.0036 

0 00*^1 + 0 00015 
U.UUOO±_ 0> 00014 

\Vcd\ 

0.224 

0.22520 ±0.00065 

\Vcs\ 

0.9736 

0.97344 ± 0.00016 

\Vcb\ 

0.0433 

0.0412±°;°oos 

| Vtd\ 

0.00853 

0.00867l“? 

\Vts\ 

0.0426 

0 0404 + °' 0011 
U.U4tU4 _ 0 0005 

\v tb \ 

0.999057 

n QQQ14B + 0 000021 
u.yyyi^D o .000046 

J 

2.98 x 10 -5 

(2.96to ig) x 10~ 5 

5 

61° 

68° 


Table III: Model and experimental values of the quark masses and CKM parameters. 

The obtained and experimental values of the quark masses, CKM matrix elements, Jarlskog invariant J and CP 
violating phase S are reported in Table IIIII We use the experimental values of the quark masses at the Mz scale, 
from Ref. ^3 (which are similar to those in fUJ), whereas we use the experimental values of the CKM parameters 
from Ref. [6|. The obtained values of the quark masses and CKM parameters are in excellent agreement with the 
experimental data, as indicated by Tabic IIIII 


V. CONCLUSIONS 


In this paper we present an extension of the minimal 331 model with /3 = — -^=, based on the extended SU(3)c <8> 

SU(2 >)l®U(1)x®Tt®Z2®Z^®Zi4 : symmetry. Our economical fa flavor 331 model, which at low energies reduces to 
the minimal 331 model with j3 = — is compatible with the experimental data on fermion masses and mixing. The 
model has in total 16 effective free parameters, which are fitted to reproduce the experimental values of the 18 physical 
observables in the quark and lepton sectors. The fa and Zj, symmetries reduce the number of model parameters. 
In particular, the symmetry determines the allowed entries of the neutrino mass matrix and decouples the SM 
quarks from the exotic quarks. The Z 2 symmetry generates the hierarchy between SM up and SM down type quark 
masses. We assumed that the SU{3)l scalar singlets having Yukawa interactions with the right handed Majorana 
neutrinos acquire VEVs at very high scale, then providing very large masses to these Majorana neutrinos, and thus 
giving rise to a double seesaw mechanism of active neutrino masses. Consequently, the neutrino spectrum includes 
very light active neutrinos as well as heavy and very heavy sterile neutrinos. We find that the heavy and very heavy 
sterile neutrinos have masses at the ~ MeV and ~ TeV scales, respectively. Thus, the MeV scale sterile neutrinos of 
our model correspond to dark matter candidates. The smallness of the active neutrino masses is attributed to their 
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scaling with inverse powers of the high energy cutoff A as well as well as by their quadratic dependence on the very 
small VEV of the Z 2 <8 <g> Z\ 4 neutral, SU(3)l singlet and T-j antitriplet scalar field The observed hierarchy of 

charged fermion masses and quark mixing matrix elements arises from the breaking of the Z 2 < 8 > Z 3 ® Z \ 4 discrete group 
at a very high energy. The tau, muon and electron masses arise from effective seven, nine and twelve dimensional 
Yukawa operators, respectively. We find for the scale of these operators the estimate A ~ 10 5 TeV. The complex phase 
responsible for CP violation in the quark sector has been assumed to come from up-type quark Yukawa terms. The 
model predicts an effective Majorana neutrino mass, relevant for neutrinoless double beta decay, with values mpp = 
3 and 40 meV, for the normal and the inverted neutrino spectrum, respectively. In the latter case our prediction 
is within the declared reach of the next generation bolometric CUORE experiment (4?} or, more realistically, of the 
next-to-next generation tone-scale O^/3/3-decay experiments. Furthermore, a vanishing leptonic Dirac CP violating 
phase is predicted in our model. 
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Appendix A: The product rules for T7 


The group X 7 , which is a subgroup of SU(3) and A(3N 2 ) with N = 7, has 21 elements, is isomorphic to Z 7 xi Zj, 
and contains five irreducible representations, i.e., one triplet 3, one antitriplet 3 and three singlets lo, li and 1 2 E3- 
The discrete group X 7 is the minimal non-Abelian discrete group having a complex triplet. The triplet and antitriplet 
irreducible representations are defined as (l 0 | : 



^1 \ 

1 X-X \ 

( x 6 \ 

3 = 

X 2 1 > 3 = 

X-2 = 

X 5 


w 

\*-4 / 

U/ 


The product rules for triplet and antitriplet tensor irreducible representations are given by: 




i ^2 (xiy6 + u k x 2 y5+u 2k x 4 y3)i k - 

fc—0,1,2 


Whereas the tensor products between singlets are: 


{x)i 0 (y)i 0 
{x)i 2 {y)i 2 


(x)i i(y)i2 = (x)i 2 (y)h = (xy) 1 0 , 
(xy) i 2 , 

(xy) ii- 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 
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The product rules between triplets and singlets satisfy the relations: 


( 2 /) i fc ® 

^ 2T(6) ^ 

x 2(5) 

- 

( 2/24(6) \ 
U x 2{5) 


\ *4(3) ) 

3(3) 

V 2/274(3) J 


(A6) 


where to = e 1 ^ . The representation lo is trivial, while the non-trivial li and I 2 are complex conjugate to each other. 
Some reviews of discrete symmetries in particle physics are found in Refs. [a-flOl l 54 j. 
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